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ABSTRACT 

We consider quasi-stationary two-dimensional magnetic reconnection in a partially ionized 
incompressible plasma. We find that when the plasma is weakly ionized and the collisions between 
the ions and the neutral particles are significant, the transition to fast coUisionless reconnection 
due to the Hall effect in the generalized Ohm's law is expected to occur at much lower values 
of the Lundquist number, as compared to a fully ionized plasma case. We estimate that these 
conditions for fast reconnection are satisfied in molecular clouds and in protostellar disks. 

Subject headings: magnetic fields — reconnection - molecular clouds - protostellar disks 



1. Introduction 

Magnetic reconnection plays a very important 
role in astrophysical plasmas. During the recon- 
nection process magnetic energy is converted into 
plasma kinetic energy, thermal heat and acceler- 
ation of charged particles, and t he topology of 
magnetic field lines is rearrang ed (lKulsrudll200l 
Yamada et al.|[2010[ ) . Magnetic reconnection is be- 
lieved to be the power source behind various as- 
trophysical phenomena, such as solar flares and 
geomagnetic storms. Magnetic reconnection also 
frequently controls transport of charged particles 
and heat in interstellar and intergalact ic media 
(iKulsrudlbodsl: IZweibel fc Yamadall20Q9h . 

In order for reconnection to be the energy re- 
lease mechanism in transient p henomena such 
as solar flares, it must be fast (iKulsrud l2005t 



Uzdenskvl l2007t lYamada etall |2010|) Although 



slow reconnection is well explained by the Sweet- 
Parker model for r econnection i n highly conduc- 
tive, hot plasmas ( Sweetl Il958l : IParkeil [l963[) . a 
common theoretical picture of fast magnetic re- 
connection has not emerged yet. A possible reason 
is that physical processes able to enhance dissipa- 
tion in a reconnection layer and to cause fast re- 
connection are fairly complicated for a theoretical 
or experimental study. However, with develop- 



ment of supercomputers, considerable progress in 
understanding possible mechanisms of fast mag- 
netic reconnection has been achieved by mea ns of 
numerical simulations ( Yamada et al. 20101) . In 
particular, one of the most important results that 
has been found both in simulations and in labora- 
tory experiments is that in fully-ionized plasmas 
the transition from slow to fast reconnection oc- 
curs when the Sweet-Parker reconnection layer 
thickness becomes comparable to the ion inertial 
length, so that the Hall term in the generalized 
Ohm's law becomes import ant (for examp l e, see 
Ma fc Bhattacharieel [19961: iBiskamp et all Il997t 



2006 



Birn et al.l200lHC assak e t al.ll2005l:lDrake fc Shav 



Yamada et a l. 2006). Since this condition is 



equivalent to the collisional mean free path exceed- 
ing the lengt h of the current sheet mu ltiplied by 
(me/m,)^/^ (|Zweibel fc Yamadal [20091 ) ■ this type 
of reconnection is also referred to as fast collision- 
less reconnection. CoUisionless reconnection can 
occur in space plasmas, in the solar corona, and in 
hot accretion disks. It cannot occur in the inter- 
stellar medium, however, unless the reconnection 
layer is very short co i npare d to macroscopic scales 
(jZweibel fc Yamadal (|2009l ): see Section i]). 

In many astrophysical systems, such as much 
of the interstellar medium and the solar chro- 
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mosphere, the ionization fraction is low. Stud- 
ies of collisional reconnection in partially ion- 
ized gases have revealed two effects which bear 
on the reconnection process. If ion- neutral col- 
lisions are sufficiently weak, the plasma and the 
neutrals decouple, so the reconnection speed is 
scaled by the Alfven speed Vai = B/^/Anpi in 
the plasma alone, while in the strongly colli- 
sional case the relevan t speed is the bulk Alfven 
speed Va = B/y/Mp (|Zweibel|[l989l ). Thus, re- 
connection with weak friction is faster than re- 
connection with strong friction by p/ Pi- A 
separate ef fect is the thinning of magne tic neu- 
tral sheets (jBrandenburg fc ZweibellHoM ). which 
can dramatically increase their merging rate 



( Heitsch fc Zweibe]|2003t iLazarian. Vishniac. fc: "Ch3 



20041: iHillier. Shibata. fc Isobell2010l) . It goes (al- 
most) without saying that neutrals also affect the 
reconnection process by making the plasma more 
resistive. This is an important effect in the low 
chromosphere, in protostellar disks, and in the 
densest interstellar gas. 

The onset of collisionless, or Hall, reconnection 
in partially ionized gases has not yet been exam- 
ined. In this paper we derive the condition for a 
transition to fast reconnection in partially ionized 
plasmas and apply our results to reconnection in 
molecular clouds, protostellar disks, and the solar 
chromosphere. We anticipate that flares in proto- 
stellar disks will be observable with ALMA, and 
chromospheric flares with IRIS, making our results 
especially timely. 

Our main results are as follows. When the Hall 
effect can be neglected, there are three regimes of 
reconnection in a weakly ionized gas, which we re- 
fer to as weak, strong, and intermediate coupling, 
respectively. When the ion-neutral collision fre- 
quency Vim reconnection current layer length L, 
and plasma Alfven speed Vai satisfy the inequality 
VinL/VAi < 1, the neutrals are decoupled from the 
reconnection process, and the reconnection rate 
is determined by the plasma parameters. When 
the neutral-ion collision frequency i>ni = t^inPi/p, 
layer length L, and bulk Alfven speed Va sat- 
isfy VniL/VA > 1, the neutrals move with the 
plasma and the reconnection rate is determined 
by the bulk parameters. In the third regime, 
the coUisionality is intermediate, and dissipation 
by ion-ncutral friction is especia hy stron g . Th e 
first two cases were described in IZweibell (|l989f ). 



and all three are analogous to the regimes of 
MHD wave propagation in partially ionized gases 
(iKulsrud fc Pearcyil969f ). 

The onset of Hall reconnection occurs in the 
weakly coupled regime just as it would for a 
fully ionized plasma (with resistivity modified 
by electron- neutral collisions). In the strongly 
coupled regime, however, while the Sweet-Parker 
layer is thickened by enhanced resistivity and 
reduced effective Alfven speed, the ion inertial 
scale can be increased even more, enlarging the 
parameter space for fast, Hall mediated recon- 
nection. A similar enhancement of the Hall ef- 
fect in weakly ionized systems has been seen 
in studies of the mag neto-rotational instability 
in protostellar disks iBalbus fc Terguem 2001 
Salmeron fc Wardlell2005[l" 



In the next section we present basic three- 
fluid magnetohydrodynamics (MHD) equations 
for partially ionized plasmas. In Section |3] we dis- 
cuss physical conditions in the interstellar medium 
(ISM). In Section m we derive equations that de- 
scribe quasi-stationary magnetic reconnection in 
partially ionized plasmas. In Section[5]we find the 
solution of the equations and analyze it. Finally, 
in Section [6] we apply our results to magnetic re- 
connection in molecular clouds, protostellar disks, 
and the solar chromosphere. 

As in many other studies of reconnection, we 
concentrate on the 2D case. We believe the three 
regimes of MHD reconnection discussed here - cor- 
responding to weak, strong, and intermediate ion- 
neutral coupling - are robust, as these have been 
found to describe many other MHD phenomena 
in weakly ionized media. The increase in the ion 
- electron decoupling scale due to increased effec- 
tive ion mass is probably similarly robust. How- 
ever, there are undoubtedly effects intrinsic to 3D 
which all 2D studies miss, and our work shares 
these limitations. 

2. Basic three-fluid MHD equations 

In this paper, except in the next section, we use 
the Hcaviside-Lorentz rationalized physical units, 
in which the speed of light c and four times tt 
are replaced by unity, i.e. c — >■ 1 and 47r — >■ 1. 
In order to convert our equations to the Gaus- 
sian centimeter-gram-second (CGS) units, the fol- 
lowing substitutions should be made: magnetic 
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field B B/V47r, electric field E cE/V^tt, 
electric current j — )■ ^/Anj/c, electrical resistiv- 
ity 77 — )■ 77C^/47r, and the proton electric charge 
e — > \/4n e/c. 

Let us consider a partially ionized, non-relativistic, 
quasi-neutral, incompressible three-component 
plasma, which is composed of electrons, single- 
charged ions, and neutral particles. The mo- 
men tum equations for these three components 
are (Braginskii 1965t Sturrock 1994 ) 



= -VPe - ne{E + u'' x B) 

- PeVeiiu" - U') - PeVeniW - U^), (1) 
[dtU' + (u'V)u*] = -VPi 
-I- ne(E -I- X B) 

+ PeVetiu" - U') - p^V,n{u' - u"), (2) 

Pn [9tu" + (u"V)u"] = -VP„ 

+ P,Z/„(U* - U") + Pel^eniu" - u"). (3) 

Here, for simplicity, we neglect electron inertia on 
the left- hand- side of eq. ([U ; n is the electron num- 
ber density, equal to that of the ions in a quasi- 
neutral plasma; u, p and P are velocity, mass den- 
sity and pressure respectively; we assume the pres- 
sure tensors are well approximated by scalars for 
all species. The subscripts and superscripts "e", 
"i" and "n" refer to the electrons, ions and neu- 
tral particles. The last two terms on the right- 
hand-side of each of the eqs. (II])-(|3|) represent the 
momentum exchange between the plasma compo- 
nents due to electron-ion, electron-neutral and ion- 
neutral collisions with effective frequencies i^ei, i^en 
and i^in respectively. For simplicity, we neglect 
electron-electron, ion-ion and neutral-neutral col- 
lisions and the corresponding viscous forces. Also, 
in this study we neglect ionization and recombina- 
tion processes, and, therefore, the densities 



nrrie 



nnii 



(4) 



of the three plasma components are constant in 
the incompressible plasma case. Here me, rrii and 
rUn are the electron, ion and neutral masses re- 
spectively, and n„ is the neutral number density. 

The electric current is j = ne(u' — u'^), and, 
therefore, the electron velocity is 

u'' = u* - j/ne. (5) 

Substituting eq. ([5]) into eq. ([T]), we obtain Ohm's 



law 

E = ?7j - X B + (l/ne)j X B 

- (l/ne)VFe - (meiye„/e)(u* 

where j x B/ne is the Hall term, and 
4 = (me/ne2)i/2. 



u"), (6) 



(7) 



Here r]ei is the standard Spitzer resistivity ([Sturrock 
|1995), ?7e„ is the resistivity due to the electron- 
neutral collisions, 77 is the total resistivity, and 
is the electron inertial length. The total resistiv- 
ity r] is enhanced over the Spitzer value by the 
electron-neutral collisions, as one expects. 

It will sometimes be useful to work with the 
electron collision time Te = {i^ei + Ven)~^ ■ In terms 
of r,., 



77 ; 



(8) 



Next, we take the sum of equations fT])-© and 
use formula ([S]). As a result, we obtain the mo- 
mentum equation for the ions: 

[dtu' + (u'V)u'] - -V(Pe + P,) + j X B 

- {Pil^2n + Pel^en){u' - u") + {mel^en/e)i, (9) 

Finally, we substitute eq. ^ into eq. ^ and 
rewrite the momentum equation for the neutrals 
as 

p„ [a^u" + (u"v)u"] = -vp„ 

+ {Pi>^i7i + Pe>^en)iu' - u") - (777e7/e„/e)j. (10) 

Equations (|6|)- (fT0|) together with the Maxwell 
equations are the basic three-fluid MHD equations 
for a partially ionized plasma. In addition to these 
equations, we note that in incompressible and non- 
relativistic plasmas the velocities and the electric 
current are divergence-free, V • u* = V • u" = 
and V • j = 0. 

3. Numerical Expressions for Parameters 

Let us estimate the values of physical parame- 
ters in representative astrophysical plasmas. This 
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is useful to motivate some approximations. In Sec- 
tion [6l we will apply the theoretical results of this 
study to the weakly ionized interstellar medium 
(ISM), protostellar disks, and the solar chromo- 
sphere. In this section we temporarily use the 
Gaussian centimeter-gram-second (CGS) physical 
units. 

The ion-neutral, electron-neu tral and electron- 
ion coUisional fre quencies are ( Braginskii 1965t 
Draine et aLlll983[) 



-I 1 /A — Q —1 Pn Pi 

i/i„« 1.9x10 sec - — , — n. 



1 + rrin/mi 
Pi rrin 

-1 „ Tn-3/2 



K ' 



(11) 



I', 



60 sec nT; 



K 



Here the frequencies are measured in inverse sec- 
onds (Hertz), the electron number density n is 
measured in cm~^, and the electron temperature 
T is in Kelvins. We multiply the expression for vin 
by a factor of four in making estimates for the solar 
chro mosp here, due to i ts rela tively large tempera- 
ture lOe^Pontieu et al.ll200l[ ). The total electrical 
resistivity of the magnetic field, given by eq. ([7]), 
is relatively small. 



2.4 X lO^^sec X T^^^^ 



X [l + 1.4 X 10-"(p„/pO("i^/m„) T\ 



As a result, the characteristic Lundquist number 
Si = VAiLgxt /ivc'^ /^'^) is very large in cosmic 
plasmas, 

5*, w 2 X lO'"^ (mp/TO,)i/2 



Lext,AU Bi,xt,p,G Tp/ '^crJ-'^ 



> 



1 + 1.4x 10-ii(p«M)(m,/TO„) Tl 
> 1. (12) 

Here Lext.AU is a characteristic system size in the 
astronomical units (AU), Bext,iiG is the reconnect- 
ing magnetic field in microgauss (pG), nip is the 
proton mass, and velocity V^i = B^xt I ^/iTrnmi is 
the Alfven velocity based on the ion density. 



VA^ ~ 2.2 X 10 — 1/2 ■ 

^^'^ my n 

t rm ^ 



(13) 



As we shall see below, it is useful to introduce 
the ion inertial length di = {miC^ /Anne^Y^^ [in 



Heaviside-Lorentz units, di = (rrii/ne^Y/'^; seethe 
last of eqs. Its approximate value is 

« 2.3 X lO^cm (m^/mp)i/2 n'^X- (14) 

Useful alternative expressions for Si and Vai 
are 

= -E^l^^^Te, VAi = —UJcedt, (15) 

di nii 

where tOce = eB/mgC is the electron cyclotron fre- 
quency. 

Using eqs. (O and ((TT|) . let us estimate the fol- 
lowing important dimensionless ratios, which show 
the relative strength of particle collisions: 



niiViji 
Vet 



0.00024 T, 



1/2 



K 



« 1, (16) 



^ « 1.4x10-"^— Tl. (17) 

Vei Pi rUn 



rrii 



Equation p6)) is an estimate for the ratio of the 
last two terms in eq. ([3]). We see that it is very 
small, unless the electron temperature is several 
millions degrees. Thus, due to relatively small 
electron mass, the effect that the neutral particles 
experience from their collisions with the electrons 
is typically negligible as compared to the effect 
from the ion-neutral collisions. In contrast, the 
relative strength of the electron-ion and electron- 
neutral collisions, given by eq. (jl7p . can be ei- 
ther large or small, depending on the density ratio 
Pn/ Pi and on the electron temperature T . There- 
fore, the total resistivity 77, given by eq. can be 
dominated by either electron-ion collisions or by 
electron-neutral collisions in interstellar medium 
and in laboratory plasma experiments. 

4. Reconnection equations 

In this section let us derive equations that de- 
scribe the magnetic reconnection process in par- 
tially ionized plasmas. 

It turns out that, when inequality rrieVen ^ 
rriiVin holds in a system undergoing magnetic re- 
connection [refer to eq. (dH)], the electron-neutral 
collisions can be neglected in all equations, except 
in eq. ([7]) for the total resistivity. The proof is 
given in Appendix [Xj As a result, we can omit 
the terms proportional to the electron-neutral col- 
lision frequency Vf-n in eqs. (jS]), (0) and (ITUl) . and 
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(2A,x2Lext 



Fig. 1. — The geometry of the reconnection layer, 
with common notations. The point O is a mag- 
netic X-point. The fluid flows toward O along the 
X axis and away from O along the y axis, carrying 
the magnetic field. The field is frozen into the elec- 
tron fluid everywhere except inside the reconnec- 
tion current layer, which has width 26 and length 
2L. In MHD reconnection the ion and electron 
decoupling regions coincide, but in Hall reconnec- 
tion the ion decoupling layer is larger; its width 
and length are 2Aj and 2Lext, respectively. The 
points O, M, and M are defined in Appendix B. 
For additional explanation, see the text. 

can rewrite these equations as 

E = 77j - u* X B j X B/ne- VPe/ne, (18) 
[dtu' + (u*V)u*] -V{Pe +PO+ixB 

- p.uUxi' - u"), (19) 
p„ [dtvT + (u"V)u"] = -VP,, 

+ p,i.„(u^-u"), (20) 

where the total resistivity rj is given by eq. ([7]) 
and includes a contribution from electron-neutral 
collisions. 

Let us now describe the reconnection layer, 
shown in Figure [TJ We assume the classical two- 
dimensional Sweet-Parker-Petschek geometry for 
the reconnection layer. The layer lies in the x-y 
plane of the coordinate system, and the x- and y- 
axes are chosen to be perpendicular to and along 
the reconnection layer respectively. The z deriva- 
tives of all physical quantities are assumed to be 
zero. 

The thickness of the reconnection current layer 
is 2(5, which can be formally defined by fitting the 
Harris sheet profile {Bext/ 5)cosh~'^{x/ 5) to the 
current profile jz{x,y = 0). The length of the 
reconnection current layer is 2L. Outside the re- 
connection current layer the z-component of the 



Ohm's law p8]) reduces to = —{u'^ x B)^ [see 
also eq. ([5])], and, therefore, the magnetic field 
lines are frozen into the electron fluid. Thus, the 
reconnection current layer coincides with the elec- 
tron layer, which is the region where the electrons 
are decoupled from the field lines. 

The ion layer, which is the region where the ions 
are decoupled from the field lines, can be much 
larger. We use notations 2Ai and 2Lext for the ion 
layer thickness and length, where L^xt is also ap- 
proximately equal to the external (global) scale of 
the magnetic field. We have Ai> 6 and L^xt ^ L. 
The region where the neutral particles are decou- 
pled from the ions can be still larger than the ion 
layer. 

The value of the reconnecting field By in the 
upstream regions outside the reconnection layer 
(at X 5) is approximately equal to the value of 
the external (global) magnetic field Bext outside 
the ion layer, up to a factor of order unity. This 
can easily be seen from the definition of S and 
from the Ampere's law 2;-component By{x,y = 
0) J^jz{x',y = 0)dx' . The out-of-plane 



ture torake & Shav 


I2OO6I: lEastwood et al. 


2007 


Yamada et al. 2010 


; Zweibel & Yamada 


2OO9I). 



Finally, the reconnection layer is assumed to have 
a point symmetry with respect to its geometric 
center, point O shown of Figure [TJ As a result of 
reflection symmetries with respect to the x- and 
y-axes, the x-, y- and 2;-components of u, B and 
j have the following symmetries: Ux{ix,^y) = 
±Ux{x, y), Uy{±x, =Fy) = TUy{x, y), Uz{±x, Ty) = 
Uz{x,y), Bx{±x,Ty) = TBxix,y), By{±x,Ty) = 
±By{x,y), Bz{±x,Ty) = -Bz{x,y), jx{±x,Ty) = 
±jx{x,y), jy{±x,Ty) = Tjy{x,y) and j^(±x,Ty) = 
jz{x, y). Here u is the velocity of any species. We 
extensively use these symmetries in the forthcom- 
ing analytical derivations , which are similar to the 
derivations in iMalvshkinI (|2008l ). 

Let us list the assumptions that we make for the 
reconnection process in a partially ionized plasma. 
First, as we have already stated above, we neglect 
ionization and recombination processes. Second, 
we assume that the collision frequencies and resis- 
tivities 77, rjei, Ven are constant in space and time. 
We also assume that the characteristic Lundquist 
number Si is very large. 

Si ^VAiLext/V^ VAi ^ Bext/^/Pi, (21) 
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an assumption easily satisfied in cosmic plasmas 
[see eq. (fT2|) ]. Note that the Alfven velocity VAi 
is calculated by using the ion density pi and the 
reconnecting magnetic field value B^xt ■ Third, we 
assume that the reconnection process is station- 
ary or quasi-stationary, so that all time deriva- 
tives can be neglected in all equations. This as- 
sumption means that the reconnection rate is slow 
sub-Alfvenic, <C VAiB^xt, and that there are 
no plasma instabilities in the reconnection layer. 
Fourth, we assume that the reconnection layer is 
thin, S <^ L and Ai <C L. This assumption is 
related to the previous assumption of slow recon- 
nection because of the mass conservation condition 
for the plasma. 

Before we proceed with derivations of the re- 
connection rate it is convenient to introduce the 
following dimensionless parameters: 

P = Pn/Pr, (22) 

i> = /2{dyui), , (23) 

V ^ {dyU-y)o/{dyV}y)o. (24) 

7 EE {dxyBz)o /ne{dyUy)o ■ (25) 

In eq. (P^ . p is the ratio of the densities of the 
neutrals and the ions. Equation p3)) introduces 
the ion-neutral collision frequency, normalized by 
two times {dyUy)o- The latter is the ion accelera- 
tion rate dyUy calculated at the central point O 
(see Figure [T]). The parameter v, defined by 
eq. p4|) . is approximately the ratio of the acceler- 
ation rates of the outflowing neutrals and the ions 
inside the reconnection current layer. Finally, the 
parameter 7, defined by eq. (j25p . is the normalized 
value of {dxyBz)o that is the second order mixed 
derivative of the quadrupole out-of-plane field B^ 
at the central point O. This parameter 7 gives 
the approximate ratio of the Hall term j x B/ne 
and the — u' x B term that enter Ohm's law (|18p . 
inside the reconnection current layer. Thus, the 
Hall term is important when 7 > 1. Note that all 
parameters (I22 p - (l25|) are non-negative. We will of- 
ten replace derivatives by inverse length scales in 
estimating these parameters. 

Now let us use eqs. p^ - (l20)) . the Maxwell equa- 
tions and the incompressibility relations 

dxv^ - ~dyv}y, dxK = (26) 
to derive the formulas that we will later solve for 



the reconnection rate and other physical quanti- 
ties. 

First, we use Ampere's law. The displace- 
ment current can be neglected in a non-relativistic 
plasma, therefore, we find 

3x dyBz, jy dxBz, ('97'\ 
i — B B - B B ^ ' 

We can estimate the ^-component of the electric 
current, jz, at the central point O as 

jo = {jz)o^ {dxBy~ dyBx)o~ 

~ {BxBy)o~B,xt/S, (28) 

where we use {dyBx)o ^ {dxBy)o ~ B^xt/S at the 
point O. The last estimate, {BxBy)o ~ B^xt/S, 
follows from the fact that S is defined as the half- 
thickness of the jz profile across the reconnection 
layer. 

Faraday's law V x E = —dtB for the x- and y- 
components of the magnetic field gives dyEz = 
—BfBx = and dxEz = BtBy ~ 0, where 
the time derivatives are neglected in the case of 
quasi-stationary reconnection. As a result, the z- 
component of the electric field Ez is constant in 
space, 

Ez = riiz - « - jx/ne)By + (uj, - jy/ne)Bx 
= constant. (29) 

Here we use Ohm's law to find the expression 
for Ez- The reconnection rate is the rate of de- 
struction of the magnetic fiux, — Bt By dx — 
- dtBy dx ^ ~ dxEz dx = {Ez)o- Thus, it 
is given by the value of Ez at the central point O, 

Ez = vjo. (30) 

One of our goals is to find the value of the recon- 
nection current jo and to calculate the reconnec- 
tion rate given by eq. pO]) . 

With the time derivatives neglected, the z- 
components of the momentum equations (|19|) 
and ^ are 

P,{ulBxUl + ulByUl) = jxBy-jyBx 

~PiVin{ul-u'l), (31) 
PniUxdxU" + UyByU") = PiVin{u\- vJl). (32) 

Calculating jdx^ of these equations at the cen- 
tral point O, dividing the resulting expressions by 
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2pi{dyUy)o, and using eqs. ([22|) - (l27t . we obtain 

(1 - l'){dxxUl)o + vidxxuDo = 

= ~j{e/mi){dxBy)o, (33) 
HdxxuDo + {~pv ~ v){dxxK)o = 0. (34) 

Similarly, taking jdy^ of equations and (15^ 
at the point O, we find 

^ j{e/mi){dyBx)o, (35) 
- HdyyuDo + {pv + v){dyyU^:)o = 0. (36) 

Equations p3p -(|36 |) are two systems of two linear 
equations in each system for the unknown quan- 
tities {dxxuDo, {dxxK)o, {dyyuDo and {dyyU^)o- 
The solution is 

{dxxOo - -7(1 + v~pv / D^){elm,){dxBy)oX'i7) 
{dxxOo = i{vlD^){elmi){dxBy)o, (38) 
{dyyuDo = 7(1 - vpv/ D+){e/mi)(dyBx)o, (39) 
(a^^Oo = ^{VlD+){e/mi){dyBx)o: (40) 
i:>± = pv±v{l + pv). (41) 

The determinant -D+ of system (|35l) - (|36|) is always 
positive, except for the trivial case when there are 
no collisions with the neutral particles [v = 0). 
At the same time, the determinant D_ of sys- 
tem ([33|) -(p4 | can be zero even when v ^ Q. For 
now we will assume that D- is non-zero, the op- 
posite case will be discussed belowQ 

Next, we calculate the second-order derivatives 
jdx^ and /dy^ of equation ([29]) at the central 
point O, using the fact that is constant. We 
find 

= ■n{dxxjz)o~'2[{dxu''^)o- (dxjx)o/ne]{dxBy)o. 

= ■n{dyyj;,)o + 2[{dyUy)o- {dyjy)olne]{dyBx)o- 

We rewrite these formulas by using eqs. (O, (j^ - 
p7| and dSHl-dlOl). We obtain 

-Tl{dxx].)o - 2{dyv'y)o{dxBy)o{l+i), (42) 

-ll{dyyJ,)o = 2{dyV'yUdyBx)oil+l). (43) 

Taking the ratio of these two equations, we find 

{dyBx)o = {^xBy)^^^ - %^ - if ' (44) 



^ We shall see that the assumption Z)_ 7^ is satisfied if 
!> <g; 1 or if i> 2> 1, the case /> 1 will be considered 
separately. 



where we use eq. (l28|l and the estimates {dyyjz)o ~ 
-jo/L"^ and (dxxjz)o ~ -jo/^'^- We again use 
these estimates and eq. (pS)) . to rewrite equa- 
tion ([H]) as 

rijl^Bl^,(dyV^yUl + ^). (45) 

Here and below we neglect all factors of order 
unity. It is noteworthy that equation (I45p de- 
scribes the supply of magnetic energy B^^i into 
the reconnection layer, where it is dissipated by 
the Joule heating rjj'^ . The rate of magnetic energy 
supply, {dyul)o{l + 7) = {dyul)o - idyjy)o/ne = 
{dyUy)o, is equal to the electron velocity deriva- 
tive because magnetic field lines are frozen into 
the electron fluid outside the reconnection current 
layer. 

Now we use Faraday's law V x E = — StB for 
the z-component of the magnetic field. We have 
dxEy — dyEx = —dtBz = 0, where the time deriva- 
tive is neglected again. We substitute Ex and Ey 
into this formula from Ohm's law (|18p and obtain 

= vidxjy - dyjx) + {Bxdxjz + Bydyjz)/ne 
+ uldxBz + UydyBz - Bxdxui - Bydyui. 

Calculating the d^/dxdy derivative of this equa- 
tion at the central point O and using equa- 
tions ([26]), ([271), dSZ]) and jSHl), we obtain 

= -■n[{dxx+ dyy){dxyBz)]^ 

+ [{dxxjz)o(dyBx)o + {dyyjz)o{dxBy)o]/ne 

+ {2ji>p^dyD^D+){e/mi){dxBy)o{dyBx)o 
« rijne{dyul)o/S^ 

-[{jo/S'){dyBx)o + ifjL')]/ne 
+{ji>p^i,yD^D+){e/mi)j,idyBx)o- (46) 

To derive the final approximate expression, we use 
the estimates {dyy)o ~ — ^ {dxx)o ~ 
and {dxBy)o ~ jo, we also use equations (|22|) - 
(|25p and we drop factors of order unity. In Ap- 
pendix |A] we show that the last term in eq. (|46|) 
can be neglected (assuming £)„ ^ 0). As a result, 
dropping this term and using eqs. (|2T|) . ([28]) . (04]), 
we obtain 

- _ _^ ^47^ 

rin'^e'^L'^{dyUy)o 

Next, we consider the acceleration of the 
plasma in the ^/-direction, along the reconnection 
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layer. We calculate d/dy of the y-components of 
the momentum equations (fTO]) and ((20|) at the cen- 
tral point O, and, neglecting the time derivatives 
for a quasi-stationary reconnection, we obtain 

P^{^vV^v)l = -{dyy[Pe+P^)o+jo{dyB^), 

- PiVin{dyUy~dyUy)o, (48) 

+ Pii'in{dyul-dyu'^)o, (49) 

In Appendix iBlwe estimate the pressure terms and 
find that 

+ 0{p^i^yUl)l} + 0{p,iy,nidyUl-dyU';)o},m 
{dyyPn)o = 0{pn{dyU'^)l} 

+ o{piiyin{dyul-dyUy)o}, (51) 

where o{...} denotes terms that are small com- 
pared to the expression inside the brackets {...} 
in the case of a thin reconnection layer {S ^ L 
and Ai < L). We substitute eqs. ([^ and (|5T|) 
into eqs. p8)) and (|49)) and then use formulas ((22)) - 
(PS)) . (P?]) . As a result, equation (|^ becomes 

pu^ = 2£>(1 - {}), (52) 

while the sum of eqs. (|48l) and ([49]) gives 

(a,«;)2(l + pw2)^^^y^2^ (53) 

where we use eqs. (PT|) and (gl]). Equation ([55)1 
describes the increase of the total kinetic energy 
of the ions and neutrals due to the work produced 
by the pressure and magnetic forces during the 
plasma acceleration in the downstream regions. 
Note that parameter v, given by eq. ([M]) . must be 
non-negative (to be more precise, < u < 1 must 
hold) because the neutral particles are dragged by 
collisions with the ions [see eq. ^^U^]. Therefore, 
the physically correct solution for v of quadratic 
equation ([5^ is 

V = { \J + 2/5i> — v) Ip w min { \/ v/p, 1 } , (54) 

where the final expression is a convenient simple 
estimate for v. 

Let us now estimate the thickness of the ion 
layer (see Figure [T]) . Note that in the upstream 
region just outside the ion layer, at a; ~ and 



y = 0, the electrons and ions are coupled together, 
the electric current is weak, the magnetic field lines 
are frozen into the electron-ion fluid, and eq. (1^^ 
reduces to — —u\.By « —u\.Bf.xt, where u\. « 
{dxuDoAi = - {dyul)oAi. Thus, we have 

A, « E,/{dyu'y)oB,,t, (55) 
Vr ~ {dyu'y)oA,^ E,/Bexu (56) 

where Vr = |u^| « {dyUy)oAi is the reconnection 
velocity. It is the velocity with which magnetic 
field lines and magnetic energy are carried by the 
plasma into the reconnection region. 

Next, let us consider the thickness A„ of the 
region where the neutral particles are decoupled 
from the ions. If the neutrals and ions strongly 
collide and move together, v = 1, they are cou- 
pled everywhere, and A„ is not defined. If the 
neutrals and ions are not fully coupled and w < 1, 
a reasonable definition of A„ is based on the lo- 
cation upstream where the inflow velocities of the 
neutrals and ions become comparable. Namely, 
ss m!^ at a; w A„ and y = 0. Unfortunately, 
we cannot estimate A„ defined this way by us- 
ing our local, analytical approach. This is because 
the profile of the neutrals inflow velocity it" as a 
function of x is unknown in the upstream region 
outside the ion layer (i.e. at x S [A^, A„])ll Find- 
ing this profile requires full solution of the govern- 
ing PDEs. Instead, we suggest a simple estimate 
A„, based on a dimensional analysis, as follows. 
First, let us note that the effective collision fre- 
quency for the neutrals is Vni = Vin/p, which is ob- 
tained by comparing the first and the last terms in 
eq. (|20|) . Second, the neutrals achieve their max- 
imal inflow velocity around the edge of the ion 
layer, |<(A„0)| « |(9.<)o|A, = (9,m^)oA, = 
v{dyUl)oA^ [see eqs. ^ and Now, we 

can make an estimate A„ w |u"(Ai, 0)|/^'„i w 
{pv/v)Ai w pjv M Ai/v > Ai in case w < 1 
[see eqs. (|23|) and ([M)) ]. Fortunately, the exact 
value of A„ does not directly influence the recon- 



^ Note that if the ions and neutrals are weakly coupled and if 
the ion pressure force can be neg lected, then th e ion veloc- 
ity ui oc a;-i/3 at a: 6 [A;, A„] llHeitsch fc Zw cibcl 200:^). 
However, the neutrals pressure force cannot be neglected 
because otherwise the inflow velocity of the neutrals would 
be much larger than eq. II54I I implies (to see this, integrate 
the x-component of eq. I|20| l over x G [A^, A„] at j/ = 0, 
and use A„ 2> A^). 



8 



nection rate and other important physical param- 
eters, calculated below. 

In the end of this section let us estimate the 
energy dissipation rate due to the ion-neutral col- 
lisions, which heat the ions and the neutrals. The 
dissipation rate per unit time, per unit volume 
is qin = PiVin('Vi' — u")^. H Therefore, the total 
dissipation (per unit time, per unit length in the 
z-direction) inside the upper right quarter of the 
ion layer is Qin = Jq Jq ^ (lindxdy. The flux of 
the (electro)magnetic energy supplied into the ion 
layer is given by the x-component of the Poynting 
vector, (E x B)^;. Therefore, the total magnetic 
energy supplied per unit time, per unit length in 
the z-direction, is Em ~ L\Ei x Bj^; w LE^Bext- 
The ratio of the dissipated and supplied energy 
rates is 



Qi' 



PiVv, 



LEzBext 
LEzBext 



(U^ 



,n\2 



dx dy 



K - u'ifdx 



(ul - u^fdy 



PiVr, 



v{l - vf / [1 + pv'^) 



(57) 



Here, to obtain the penultimate expression, we 



use estimates (m^ - 



n\2 



and 



[dxU 
[ul - u. 



dxK 



l2 ^2 



ri\2 



{dyU 



dyK)W = {dyul)l{l-vfy^ [seeeqs. (HI, (EH)]; 
to obtain the final expression, we use inequality 
Ai < i, eqs. (121]), dMl, dMl), dSH), and we ne- 
glect factors of order unity. Note that, due to 
eq. ([53]), Qin/£m £ 1, as one expects. We also 
see that Qin/£-m increases with collisionality for 
small z/, reaches a peak value that is around unity 
as collisionality increases, and then declines again 
as -D — ^ 1. 

5. Solution 

We solve the nine equations dH]), (|iil) . 
(|I5]) . P7|. dSU-dSS]) for nine unknowns: jo, E^, 
Vr,, (S, Ai, {dyU\)o, {dyBx)o, V and 7. For the 
presentation of the solution, it is convenient to ex- 
press resistivity 77 and ion density pi in terms of 



^ To derive this formula, add together eq. II19I I multiplied by 
u' and eq. I I20I I multiplied by u". 



Lundquist number Si, Alfven velocity Vai, field 
B^xt and scale Lea;*, see eq. ([U]). It is also helpful 
to express the ion charge density ne in terms of 
the ion inertial length di, 



ne = ^i/di = Bext/diVAi, 
di = (mi/ne^y^'^ . 



(58) 



^From eqs. (|4T]) and (|54]) . we find that Z)_ = 
ii 1/ <^ 1, and Z?_ = —pvi) if 3> 1. Thus, if 
£/ 9^ 1, then condition D_ ^ is satisfied, and the 
solution is 



m = 1 + n\va{v, p}, 

w [vtnL/VAi) [1 + Tam{viriL/VAi, 
V ~ min| \/ v / p, l} , 
7 « Sidf / LLext, 



{dyUy)o ~ VAi/L\ffii, 



Jo 



E, 



^fSiB^xt 
%/ LLext 

VAiBexf 



1 

rh 

V Lext 



{dyB, 



ext 



LLext 
1 

TO 
1 

fh 



1/2 



LLe 



1/2 



VLLe 



A,« 

^out 



%L^/-^ 

1 ^ 

TO LLext 
1 



S^dl 



LLext 
-1/2 



LLe 



V LLext r 
Ez/Bext « (A./L)<„„ 

%i {dyUy)oL w VAi/V^. 



S^d} 



1/2 



(59) 

(60) 
(61) 
(62) 
(63) 

(64) 
(65) 
(66) 
(67) 

(68) 

(69) 
(70) 



Here we define to, which is the factor by which 
the ion particle mass is effectively increased due to 
the ion- neutral collisions (as we shall see below). 
We also introduce Uq^j, which is the ion outflow 
velocity in the downstream region outside the re- 
connection layer (i.e. at a; = and y ~ L). Note 
that eq. ((69)) essentially represents the mass con- 
servation law for the ions. Let analyze the above 
solution for the case 



p ^ pn/pi > 1, 



(71) 
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which holds for molecular clouds, protostellar 
disks, and the solar chromosphere 

Depending on the value of parameter v, we have 
the following cases for magnetic reconnection. 

The first case is when £> ^ 1, and, as a re- 
sult, m = 1 in eqs. ([59 |) -(f70 |) . In this case the 
ion-neutral collisions are negligible because their 
frequency is very small compared to the ion in- 
flow and outflow rates, Vin ^ (,dyUy)o ~ Vaz/L « 
Uy{L)/L « \ul.{Ai)\/ Ai. The ion-neutral cou- 
pling is weak. The neutral particles carry a neg- 
ligible fraction of the total plasma kinetic energy, 
pv"^ w i> ^ 1. Energy dissipation due to the ion- 
neutral collisions is very small, Qm/^m <C 1 in 
eq. (|57| . 

The second case is when !> ^ 1, and, there- 
fore, m ^ 1 [assuming eq. (I7T|) holds]. In this 
case the ion-neutral collisions significantly influ- 
ence the reconnection process because their fre- 
quency is large, i/^ > Vai/L > Uy{L)/L w 
|u^(Ai)|/Ai. The neutral particles carry most of 
the plasma kinetic energy in this case, j5v^ 3> 1. 
There is significant energy dissipation Qm ~ £m 
due to ion-neutral collisions if 1 ^ < p and 
the ion-neutral coupling is intermediate {v < 1). 
However, if ^ p and the ion- neutral coupling 
is strong (0 = 1), this dissipation is negligible, 
Qin ~ {p/i')£m < £m, [scc cqs. and ([571)]. 

In the case i/ 3> 1, the ion- neutral collisions 
result in an effective increase in the mass rrii of 
the ion particles by factor m. This is because 
equations ([60]), (I62|)-(l68l), ^ can be obtained 
from the corresponding equations in which m is re- 
placed by unity, by making the following substitu- 
tions: mi -J> rhrrii, pi rhpi, Vai = Bext/^/Pi 
V ai/V^, Si = VAiLext/i] Si/\ffii and = 
—7' y/Adi. 



In the limiting case of strong coupling when ion- 
neutral collisions are extremely frequent, i) ^ p 
and rh = p, the neutral particles are well coupled 
to the ions and move together, v = 1 and u" = u* 
(also Qin <C £m)- In this case the neutrals and 
ions behave as a single fluid of density pi + pn — 
{l+p)pi ~ ppi, and {dyul)o ~ Va/L = VAi/L^ffii. 
These theoretical results are in good agreement 



^ The case p < 1 is not very interesting because in this case 
m Si 1, the effective ion mass is comparable to mi, and 
the neutrals do not significantly influence the reconnection 
process, see eqs. I|59| l- I|70|l . 



with recent numerical s imulations of reconne ction 
in solar chromosphere ( Smith fc Sa kai"2008l ) , and 
with previou s theoretic al studies (jZwcibcl 1989t 
Zaqarashvili et al. 2011 ). 

The last case for magnetic reconnection left to 
consider is when i/ w 1. In this case the determi- 
nant _D_ of the system of equations (|33|) - (p4| is 
close to zero, and we find from this system that 
{dxyBz)o C!c 7 is also close to zero. As a result, 
higher order Taylor expansion terms have to be in- 
cluded into our derivations in order to estimate the 
physical quantities inside the reconnection layer in 
a mathematically rigorous way. Fortunately, we do 
not need to go through these tedious calculations. 
Instead, we note that solution (|59p - ([7n|) is contin- 
uous a,t V 1. Therefore, the case w 1 is not 
special, eqs. (|59 l) -([70 l) still hold, and the approxi- 
mate solution in this case is similar to that in the 
weak coupling case v <^ 1 because in both cases 
TO « 1. 

Next, depending on the value of the Lundquist 
number 5^, there are two distinct reconnection 
regimes that the solution (|59|) -(|70 |) describes. 

First, when Si <^ Ll^^/dfy^ (i.e. when S 
computed from the classical Sweet-Parker the- 
ory is larger than diifi^^'^), a modified Sweet- 
Parker reconnection regime takes place, for which 
A; « (5 w m^l^Lextl^fSi, L « Lext (bccause 
A, « 5), 7 « VrhS^dyLl^^ < 1, jo w 

^/S'iBext/Lextm^^*, Ez « VAiBext/VS'im'^^'^, 

{dyBx)o « rn^^'^BextlLextVSi, and the quadrupole 
field Bz « {dxyBz)oLS = ne{dyul)ojLS < B^xt- 
The difference between this regime and th e classi- 
cal Sw eet-Parker reconnection (Sweet I Q^stlParker 
19631) is that in the former the ion particle mass 



uii is effectively increased by the factor to due to 
ion-neutral collisions. 

Second, there is a Hall reconnection regime 
when Si « Ll^^/dfy/^. With the Lundquist 
number value 5'^ w L^^^/ d^^/m substituted in, 
equations ([Sn)) - ([7n)) give A^ sa di^ffii, 5 ~ 
Ai{L/L^xt) < Ai, L < Lext, 7 « Lext/L > 1, 
jo » BextLext/diL V^, Ez « Vki-Bext (di/L) (also 
see ICowlevI Il985 . for the fully- ionized plasma 
case), {dyBx)o ~ Bextdi/LLext, and the 
quadrupole field Bz « {dxyBz)oL5 ~ B^xt is com- 
parable to the reconnecting field B^xt- 

Unfortunately, our approach does not allow us 
to calculate the reconnection layer length L in the 
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HaU regime. However, similar to Malvshkin ( 20091 
a plausible conjecture can be made that 
the above Hall reconnection regime represents a 
transition to fast collisionless reconnection, during 
which the reconnection (electron) layer thickness 
6 decreases from dp to the electron inertial length 
de- This conjecture is based on numerical simula- 
tions, theory, laboratory and space observations 
of magnet ic reconnection in fully ionized plas- 
mas (e.g.. iBiskamp et al.lll997t IShav et al.lll998t 



Pritchetdl200lt ICassak et a ll 120051 : IWveant et al 
20051: iDaughton et al.l|26o6: Drake fc Shavl l2006t' 
Karimabadi et al. ''2007: Dr ake et al.ll2008l : lji etal 



20081 : lYamada et al.ii2010i ). Note that at the onset 



of Hall reconnection the layer length is L « Lext] 
the transition to fast collisionless reconnection is 
accompanied by shrinking of L relative to Le^t- 

When the layer thickness S reaches d^, electron 
inertia effects become important. In this study we 
omitted electron inertia, and, therefore, we can- 
not describe this fast reconnection regime, which 
we plan to consider in the future. At the present 
time, the important result for an application to as- 
trophysical systems is that the Hall term becomes 
important in the generalized Ohm's law, the on- 
set of Hall reconnection occurs, and a transition 
to fast collisionless reconnection happens when 



r2 



LextV 
VMd^ 



= v/l + min{i>,p}. (72) 

Here Ssp = Lext / V^i is the classical Sweet-Parker 
reconnection layer thickness, and v is given by 
eq. (|60p with L w L^xt for the onset of Hall re- 
connection. It is important that, if 3> 1 and 
p ^ 1, then the Lundquist number value at which 
the transition to fast reconnection occurs is much 
lower than the corresponding value for the fully 
ionized plasma case i> — 0. 

Using eq. (ITSt . eq. (|72)) can be recast as 



L 



Hall 



diLxJceTeVfh 



(73) 



for the onset of Hall reconnection. Let us now an- 
alyze the onset of Hall reconnection in terms of the 
global scale Lext- Refer to eqs. ([5]), (pD)) and ([75]) . 
First, note that as the value of Lext decreases from 
large to small, the ion- neutral coupling changes 
from strong to intermediate at Lext ~ -^eajt' = 
y/pVAi/vin (when i> ^ p), and changes to weak 



coupling at Lext ~ Ll^^'" = VAi/vin (when i> sa 1). 
Now, we consider how the onset of Hall recon- 
nection depends on whether it happens in the 
strong, intermediate, or weak coupling case. If 



> 



■q/dj = {me/mi)i 



then the Hall re- 
connection onset occurs at i^^" > ^eJjtN when 
the coupling is strong and fh — p. If Vin ^ 
77/df — {me/mi)T~^, then the onset occurs at 
-^ext'' ^ L]^l" , when the coupling is weak and 
m = 1. Finally, if ~ rj/df = {me/mi)Te^, 
then the onset of Hall reconnection happens in a 
range Lg^(™ < ^f^t" ^ -^eSN which is equivalent 
to VA^d^/^J < Lgf < y^VMd^r]. Our analysis 
does not yield a more precise criterion in this case. 

6. Discussion 

Let us apply our results to magnetic reconnec- 
tion in molecular clouds, protostellar disks, and 
the solar chromosphere. 

Using eqs. pT|) - p^ . we make the following es- 
timates 



VinLext n 1 o ~ V'Tli/mp 

^«0.13p 

51p _ 



3/2 



■.L 



ext,AU 



1 + mn/rrii 



B, 



, (74) 



2.2 X 10' 




^ext 

S;dJ 



X 1 + lAxlO-''p{m^/mn)Tj 



Here, as in Section[31 characteristic scale Lext is in 
the astronomical units (AU), magnetic field Bext is 
in microgauss (fJ-G), electron number density n is 
in cm~^ , and temperature T is in Kelvins. Param- 
eter v is calculated by substituting L — Lext into 
eq. dnni) and using eq. ((7^ (recall that L w Lext 
for the onset of Hall reconnection) . 

The second column in Table [T] lists the typical 
values of ph ysical parame t ers in molecular clouds. 



taken from iMcKee et al 



( 19931 ). and the corre- 
sponding values of di, Si, p, vjp, "qenMei and 



m. ^From the values given for vjp and 
TO we see that in molecular clouds the 



ions and neutrals are typically strongly coupled, 
and condition (|72p for the onset of fast magnetic 
reconnection can be satisfied. 

For the protostellar disks, we assume the fol- 
lowing dependence of physical parameters on the 
disk radius r, which is me asured in the astro- 
nomical units (IWardldlioOTl) . The neutral num- 



ber density is n„ « 5.8 x 10 



14 , 



-11/4 



• AU 



the 
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electron number density is n 10 ^^n„, the 

— 1/2 

temperature is T = 280 r^J the character- 
istic length Lext « /i « 0.03 AU (equal to 
the disk vertical scale h), and the magnetic field 

''^Au'^ ^ (based on a theoretical estimation 
of the angular momentum transport in the disk). 
Taking 0.1 AU <r< 100 ^J7, we obtain the val- 
ues reported in the third column of Table [T] We 
see that there is a strong ion-neutral coupling at 
small radii and intermediate ion-neutral coupling 
at large radii. Condition (f72|) for onset of fast re- 
connection can again be satisfied. 

In the solar chromosphere, on the other hand, 
the width of the Sweet-Parker layer generally far 
exceeds the ion skin depth, meaning that condi- 
tions for fast Hall-mediated reconnection are un- 
favorable (see the last column in Tabled]). This 
is due to the relatively high density and high level 
of ionization compared to the other two systems 
discussed here. 

We conclude that fast coUisionless magnetic re- 
connection may indeed be possible in partially ion- 
ized plasmas in molecular clouds and in protostel- 
lar disks. 

We thank Fausto Cattaneo, Hantao Ji, Arieh 
Konigl, Eric Lawrence, Masaaki Yamada for 
useful and stimulating discussions. This work 
was supported by the NSF Center for Magnetic 
Self-Organization in Laboratory and Astrophys- 
ical Plasmas at the Universities of Chicago and 
Wisconsin-Madison. 
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Table 1: Typical values of physical parameters 





molecular clouds 


protostcUar disks 


solar chromosphere 


ions 


HCO+, 29mp 


Me+ 24m„ 


, Tflp 


neut fells 


pTr> '2tti,<^ 


TTo 2 777,^ 


TTrj 2 777,^ 




10^-10^ 


2x lO^-Sx lO^*" 


10" -10^^ 








10^1(1 + i0"''"O 


T(K) 


3 30 


30 10-' 


GOOO 




10 100 iiG 


(ixlO^^ 


1 10-' G 


Lext [AU] 


IG^ 10** 


0.001 10 


10"* 10"^ 


d^ (AU) 


10-''-5xl0-'^ 


10-^-2x10-'' 


5xl0-"-5xl0-i^ 




3x10^^-3x10^" 


10^- 10** 


3x10^-3x10'" 


p 


2 X 10'^ - 2 X 10^ 


8x10"' 


1 10* 


v 


4x10^-4x101'^ 


3x10^-10^" 


3x10^-10" 


v/p 


2-2x10* 


0.04-10* 


lO^-lO'" 




0.004-0.04 


10* 10" 


5xl0-*-3 




0.03-103 


0.001-50 


500-10" 
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A. Equations with the electron-neutral collisions included 

In this appendix we prove that, if rueh'en ^ 'm-ii^im then the terms that are related to the electron- neutral 
collisions and are proportional to i>en can be neglected in the reconnection equations, except in eq. ([7]). While 
some of these terms are clearly small, for example Pei^en ^ Pi^^in in eqs. © and (jlOp . it is not immediately 
clear that other terms can be omitted. A rigorous proof requires deriving general equations with electron- 
neutral collisions included, and then showing that these terms are negligible. The general equations can 
be useful for the case when electron- neutral collisions are significant, as long as electron inertia can be 
neglected To save space, we do not give all details of the derivations, which are tedious but straightforward 
to carry out along exactly the same guidelines that are thoroughly described in the main text. Instead we 
give only the key formulas and results. Also, below we prove that the last term in eq. (|46p can be neglected 
if £>_ ^ 0. 

For brevity of notation, we assume that spatial derivatives are to be calculated with respect to all indices 
that are listed after the comma signs in the subscripts, e.g. y = dyU^ and B^^xy = dxyB^- 

With the electron-neutral collisions included into derivation, it is convenient to replace eq. (|23p by 

i? = D, + De, Di = iyin/'^{Uy,y)o, i^e = {me/mi)i^en /2{ul^y)o , (Al) 

i>e/i> ~ Ve/Vi < 1, Ve^/v W Vel / h < 1. (A2) 

Here, v^/vi <^ 1 follows directly from rUpVpn ^ fnii^im while Vf.'^/i'i <^ 1 can be used because we are 
interested in the conditions for a transition to fast collisionless reconnection at 7 « 1 (see Section [5]). 

With the electron-neutral collisions included, eqs. ([2T|) . (|22|) . ([24l) - (|28|) . ([55]) and ([56l) are unchanged. 
Equation (|29|) becomes 



Ez = Vjz - « - jx/ne)By + (u^ - jy/ne)Bx - (mei^e„/e)(u!. - m") = constant. 
We see that equation (pO| stays the same. Equations (l37|) - (l40l) become 

iul,xx)o = + Dpv/D^){e/mi){By^^)o~ {C'ePv/D^){l/ne){j^^a:x)o, 

«,xx)o = jiD/D_){e/mi){By^x)o + ii>e/D_){l/ne){j,^x^)o, 

i'ul,yy)o = 7(1 - vpvlD+){elmi){B,^^y)o + {v^pv / D+){l/ne){j^^yy)o, 

«aa)o = l{v/D+){e/mi){B.x^y)o- {velD+){l/ne){j^^yy)o, 

where D± is still given by eq. (pij) with v defined by eq. (jAip now. We assume that Z3_ is non-zero (the 
case D_ = is discussed in Section [5]). Equations (|42]) and (|43|) become 

- [7? + Ven^eil + pv) / D , ^^x) o = 2{v^y y) o{By o[l + 7(1 + Vepv/D^)], (A3) 

- [■n-TlenVe{l+pv)lD+]{j,^yy)o = 2{v^y y) o{B , ^y) o[l + 1 {I - VePv/D+)]. (A4) 

Taking the ratio of these two equations, we obtain a general version of eq. (|44|) : 

(R \ 1 + 7(1 VePv/D-) f] - TJenj^ejl + pv) / D+ . . 

^ 1 + 7(1 - '^ePv/D+) 7? + VerMl + Pv) / D ^ ' ^^^^ 

We can rewrite eq. (|A3I) in an approximate form that corresponds to eq. (j45p . 

[7? + 77e„t>e(l + pv) / D « i?Lt(";,Jo[l + 7(1 + yePv/D_)]. (A6) 



^ Note that electron- neutral collisions result in an effective increase in the electron mass me, similar to the increase in rrii due 
to ion- neutral collisions (see Section [SI . 



14 



Next, with the electron- neutral collisions included, eq. p5)) becomes 

= -Vlidxx + dyy){B^^^y)]^- {meiyen/e)[id^x + dyy){ul y-u'ly)]o 

+ [ijz,xx)o{Bx,y)o{i + VePv/D^) + {j z .yy) oiBy ,x) o{l ~ V ^pV / D +)] / Ut 

+{2^iyp''vyD^D+){e/m,){By^,)o{B^,y)o~ 
« [777 + 77e„(l - v)]ne{u'y^y)o/S^ - [{jo/S^){B^^y)o{l + v,pv/D^) + {fjL^){l - v.Pil/ D+)]/ne 
+{^i)p%^/D^D+){e/m,)jo{B^,y)o. (A7) 
Finally, eqs. ([5^ and (15^ become 

= 2£>(1 -v) + 2i>ej, (A8) 
P.(«;,,)o(l + P^^) « SL/i' + Jo(i?x,,)o (A9) 

respectively [also refer to eqs. (jBlip and (|B12p ]. The parameter v must be non- negative, therefore, the 
physically correct solution of eq. (|A8p is 



v = il/p)Wi^^ + ^Pii> + i>el)-i>]- (AlO) 

Next, let us use inequality i'e/'^i ^ 1 or, equivalently, Ve/i^ ^ 1 to simplify the above equations. We still 
assume that given by eq. is not close to zero, and, therefore, \D±\ fa pv + D{1 + pv). We have 

VePV _ VePV < £l ^ ^ ??e«I>e(l + Pt)) ^ Ve{l + Pv) < ^ ^p^-y^^ 

\D±\ pv + i>(l + pu) ~ i> ' ?7 |Z3±| '^^^ p-u + ;/(l + pv) ^ v 

As a result of inequalities and (lAlip . equations (|X5|) . ([Xe)) . (j^B and (|A10I) reduce to equations pit . 
(|45|). (|53|) and ((54)) respectively. Henceforth, we can use the latter instead of the former, except for calcula- 
tions oil — i) for which the more accurate equation (|A10|) should be used when 1 — -D <C 1. 

Now, the only proof left to do is to show that eq. (jA7l) reduces to eq. (|Tf)) . This proof is done as follows. 
Using eqs. (HI]), dM]) and (jAllI) . we rewrite eq. (jXTl) as 

ril{y^y,y)o + Ven{l " i))(w;,,)o " B^ v? + {^v p^ v^ / D ^ D +){BI,^VIJ jI) « 0. (A12) 

Let us estimate the ratio of the last term and the first term on the left-hand-side of eq. (jA12p . 

vp^v"^ -Sg^jFj- vp^v'^ 1 + pv'^ ^ V 1 + pu^ ^ V ^ pv^ 



\D_\D^^L'^fMy,y)o \D-\D+ 1+7 - {X + vY 1 + 7 " (1 + 1^(1 + 7) 

Here we use eqs. (|45p and ([53)) to obtain the second expression; we use an estimate \D± \ ^ pv + i>(l + pv) > 
pv{l + D) to obtain the third expression (assuming I?_ 7^ 0); and we use estimates < (1 + i>)^, 1 + 7 > 1 
and pv^ < i> [see eq. ()54l) ] to obtain the final result. Thus, the last term on the left-hand-side of eq. (|A12p 
can be neglected because it is comparable to the first term or smaller. 

Next let us prove that the second term, r7e„(l — v){Uy y)o, can be neglected in eq. ()A12p as well. First, 
the ratio of this term and the term B^^i/n^e^L^ is 

n^e^?7en(l - v){ul y)oL^ _ p.meVen (1 " ^)K,y)o-^^ _ 2£^e(l " v){'U^y,y)W v{l - v) , . ^ „^ 
^ext rrii ^ext ^Ai I + pV 

Here we use eqs. (j4]) and ([7]) to obtain the second expression; we use eqs. (|2T)) and ()Aip to obtain the third 
expression; and we use eq. ([53]) and inequality ^ v to obtain the final result. Second, the ratio of the 
second and the first terms in eq. \K12\ is 

~, , / = — < (A14) 

Vl[Uy,y)o VI 1 
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In the case i) < p{l + i>e'y)/i^, from eq. (jA10[) we find v ~ + v^'^)/ p and 1 — v ~ 1. Therefore, 
u{l — v)/{l + pv^) ~ + i> + C'ej) £ 1 in eq. (jAlSp . Thus, in this case the term ?7en(l — v){ul,y)o can be 
neglected in eq. (jAl2p because this term is small in comparison with the term Bl^^/n^e^L^. 

In case v 3> p{l + VfT^^jv, from eq. (jA10[) we obtain v k. \ and 1 — w w —VfT^jv + pjv. We substitute 
1 — u sa —Vg^^/v into eq. (IA14P and find |1 — v\/^ w v^/v ^ 1. We also substitute 1 — v ~ pjv into eq. (IA13[) 
and obtain i>(l — £;)/(l + pt';^) PS p/(l + p) < 1. Therefore, the term ryen(l — t')(u5, j^)o can again be neglected in 
eq. (IA12P because this term is small in comparison with either the term r]^{Uy y)o or the term B^^^/n^e^L^. 

As a result of the above estimates, the second and the fourth terms in eq. (jA12p can be omitted. Therefore, 
this equation and eq. (IA7p reduce to eq. (H71) . 



B. Derivation of equations ()50p and (j5ip 

As in the previous appendix, here we include electron-neutral collisions, and, to save space, we again 
assume that spatial derivatives are to be calculated with respect to all indexes listed after the comma signs 
in the subscripts, e.g. y = dyU^ and B^^^y = dxyBz- 

Our derivation of equations (l50l) and (|511) , to some degree, is similar to the Sweet-Parker arguments for the 
pressure drop along and across the reconnection layer. To be precise, we integrate pressure gradient vectors 
along the rectangular contour O — > A/ — > 7\/ — > O shown in Figure [T] and use the force balance condition 
for the plasma slowly inflowing across the reconnection layer. Note that point O is the reconnection layer 
center, and point M is in the upstream region at the edge of the ion layer at x ~ A,; and y = 0. Let y be 
the y-coordinate of points M and O. We proceed as follows. 

We consider the limit when points O and M are infinitesimally close to point O and M respectively, and, 
therefore, y — > -1-0. For infinitesimally small values of the y-coordinate, we use Taylor expansions in the y 
coordinate for the x- and y-components of the velocities, current and magnetic field. 



u"J'\x) + {yy2)u:il\x), 



].=j^\x) + {y^l2)j'^M, 
Bx = yB'^lix) + {yy6)Bi''lyy{x), 



AO) 



Uy = yUy^y\x) + {y^/6)Uy':yly{x), 
< = yUyf\x) + {y^lQ)Uy'}y^y{x), 



(0) 



(Bl) 



B,, 



Bi">{x) + {yy2)B^:>y{^). 



(0) 



Here the variables with the superscripts are calculated at y = and depend only on coordinate x. 
Neglecting time derivatives for a quasi-stationary case, we rewrite the momentum eqs. (|19p and (|20p as 



V(Pe + 



B^/2) 



-pr{u'V)u' + (BV)B - {p^V^n + PeVen){n' - u") + {nieVen/e)l (B2) 
-p„(u"V)u" + {p,V,n + PeVen){Vi' - u") - (meiye„/e)j. (B3) 



Next, we calculate the line integrals of all terms in these two equations along the contour O — Af — A/ — O 
(see Fig. [1]). In these calculations we use Taylor expansions (jBl[) . and we keep only the terms up to the 
leading, second order in y (because y — )■ -|-0). As a result, we obtain 



[V(Pe + P,: + SV2)]dl = [Fe(0,y) +P,;(0,y) +i?2(0,y)/2] - [Pe(0, 0) -f P,(0, 0)] 

= (yV2)[(Pe,,,+P.,,,)o + (P.,,)?], 
P„(0,y)-P„(0,0) = (yV2)(P„,,,)o, 



j [vp„: 



dl 



;(u*V)u*]dl = (r/2 



I'M 

x,yy ) 2 / f^x,x ^x,yy 

Jo 



(yV2)o{«,)XA'} = {fmo{{v^yj,], 



(B4) 
(B5) 

(B6) 
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[(u"V)u"]dl = o{«,)^}, by analogy with eq. (Hi, 

[(BV)B]dl 



(B7) 



M 



(u^ - u")dl 



jdl 



L Jo 

{f/2) 0{iul, - <JoA./L} = if/2) o{i 

I'M 

'- Jo ^ 

{f/2)0{{B,..y)oS/L} = if/2)o{{jy,y)o} 



(B8) 



(B9) 



(BIO) 



Here, the variables inside parentheses (...)o are evaluated at the central point O; the variables with the 
superscript are calculated at point M; notation 0{...} denotes terms that, in absolute value, are 
comparable to the terms inside the brackets {...}; and notation o{...} denotes terms that are small com- 
pared to the terms inside the brackets. In derivations of eqs. (|B4|) - (|B10[) we use equations j, = —u' 



V,V' 



Jx^x 



Jv,v 



= B 



z,xy 



To derive the final approximate expressions in eqs. (|B4p - (jB10() . 
we use the following estimates at point M (which is at the ion layer edge): dy « dx ~ 



{ul^y)oA^. By 



Be 



joS, Bx 



{Bx,y)o, jz ~ Bext/L, 



jx ~ {jx,x)oS = {Bz,xy)o& (note that, in the Hall regime, the Hall term supports Ez outside of the elec- 
tron layer), and 5 < A; <C i for a slow reconnection in a thing layer. We also use the following esti- 
mates for the integral terms: J*^ Ux^x'ux^yl dx « /^«,a;)o xL~'^dx w {Uy y)lAl /L'^, {By°^ Bi°l) ^yydx w 

rM , 1,(0) ' n,{0)\ 



JoiBx,y)oBextL dx « {Bx,y)ojo6 A/ , J^\l 



)dx « So^K,x - <,x)oxL ^dx « -{ul^y 



u'^ y)oAl/L'^, a.nA Jq jx'lydx Jq Ux,x)oS L "^dx Ki -{Bz,xy)o5A,/ L'^ . 

Taking the line integrals of eqs. (|B2[) and (jB3p . using eqs. (|B4p - (jB10[) and an estimate By"^'' By"l\ « 
— i?g^t/L^, we obtain 

{Pe,yy + Pi,yy)o ~ ~BI^J + o{pi{u'yy)l] + o{jo{Bx,y)o] + o{{meVen I e){jy .y) o] 

+o{{piUiri + Pel^en){ul y -Uy y)o}, (BH) 
{Pn,yy)o = o{ p„ (u^ ,^)^ } -f o{ (mgl^en /e) (jy,y )o } + o{ (p,l/i„ -|- PeJ^en) (u^,y " ^ )o } ■ (B12) 

These equations reduce to eqs. ([50]) and ([5T|) when the terms associated with the electron- neutral collisions 
and proportional to Ven are neglected. Note that the term By™'' By"l\ « —B^^^/L^ represents the drop of 
the outside magnetic pressure along the reconnection layer. 



17 



REFERENCES 

Balbus, S.A. & Terquem, C. 2001, ApJ 552, 235 

Biskamp, D., Schwarz, E., & Drake, J. F. 1997, 
Phys. Plasmas, 4, 1002 

Birn J., Drake J. P., Shay M. A., Rogers B. N., 
Denton R. E., Hesse, M., Kuznetsova, M., Ma, 
Z. W., Bhattachargee, A., Otto, A., & Pritch- 
ett, P. L. 2001, J. Geophys. Res., 106, 3715 

Braginskii, S. I. 1965, Rev. Plas. Phys., 1, 205 

Brandenburg, A., & Zweibel, E.G. 1994, ApJ, 
427, L91 

Cassak, P. A., Shay, M. A., & Drake, J. F. 2005, 
Phys. Rev. Lett., 95, 235002 

Cowley, S. W. H. 1985, in Solar System Magnetic 
Fields, ed. E. R. Priest (Dordrecht, Holland: D. 
Reidel Publ. Co.), 121 

Daughton, W., Scudder, J., & Karimabadi, H. 
2006, Phys. Plasmas, 13, 072101 

De Pontieu, B., Martens, P. C. H., & Hudson, H. 
S. 2001, ApJ, 558, 859 

Draine, B. T., Roberge, W. G., & Dalgarno, A. 
1983, ApJ, 264, 485 

Drake, J. F., Shay, M. A., & Swisdak, M. 2008, 
Phys. Plasmas, 15, 042306 

Drake, J. F., & Shay, M. A. 2006, The fundamen- 
tals of collisionlcss reconnection, in Reconncc- 
tion of Magnetic Fields: Magnctohydrodynam- 
ics and Collisionlcss Theory and Observations, 
ed. J. Birn & E. P. Priest (Cambridge, UK: 
Cambridge Univ. Press), 87 

Eastwood, J. P., Phan, T.-D., Mozer, F. S., Shay, 
M. A., Fujimoto, M., Retino, A., Hesse, M., 
Balogh, A., Lucek E. A., & Dandouras I. 2007, 
J. Geophys. Res., 112, A06235 

Ji, H, Ren, Y., Yamada, M., Dorfman, S., 
Daughton W., & Gerhardt, S. P. 2008, Geo- 
phys. Res. Lett., 35, L13106 

Heitsch, F., & Zweibel, E. G. 2003, ApJ, 583, 229 

Hillier, A., Shibata, K., & Isobe, H. 2010, PASJ, 
62, 1231 



Karimabadi, H., Daughton W., & Scudder, J. 
2007, Geophys. Res. Lett., 34, L13104 

Kulsrud, R. M. 2005, Plasma Physics for Astro- 
physics (Princeton, NJ: Princeton Univ. Press). 

Kulsrud, R.M. & Pearce, W.P. 1969, ApJ, 156, 
445 

Lazarian, A., Vishniac, E. T., & Cho, J. 2004, 
ApJ, 603, 180 

Ma, Z. W., & Bhattacharjee, A. 1996, Geophys. 
Res. Lett., 23, 1673 

McKee, C. F., Zweibel, E. C, Goodman, A. A., 
& Heiles, C. 1993, Magnetic Fields in Star- 
Forming Regions - Theory in Protostars and 
Planets III, ed. E. H. Levy, J. I. Lunine, et al. 
(Tucson, Arizona: Univ. of Arizona Press), 327 

Malyshkin, L. M. 2008, Phys. Rev. Lett., 101, 
225001 

Malyshkin, L. M. 2009, Phys. Rev. Lett., 103, 
235004 

Malyshkin, L. M. 2010, Phys. Scripta, T142, 
014033 

Parker, E. N. 1963, ApJ Suppl. Ser., 8, 177 

Pritchett, P. L. 2001, J. Geophys. Res., 106, 3783 

Salmeron, R. & Wardle, M. 2005, MNRAS 361, 
45 

Shay, M. A., Drake, J. F., Denton, R. E., & 
Biskamp D. 1998, J. Geophys. Res., 103, 9165 

Smith, P. D., & Sakai, J. I. 2008, A&A, 486, 569 

Sturrock, P. A. 1994, Plasma Physics (Cambridge, 
UK: Cambridge Univ. Press). 

Sweet, P. A. 1958, in Electromagnetic Phenomena 
in Ionized Gases, ed. B. Lehnert (New York, 
NY: Cambridge Univ. Press), 123 

Uzdensky, D. A. 2007, Phys. Rev. Lett., 99, 

261101 

Wardle, M. 2007, Astrophys. Space Sci., 311, 35 

Wygant, J. R., et al. 2005, Geophys. Res. Lett., 
110, A09206 



18 



Yamada, M., Kulsrud, R., & Ji, H. 2010, Rev. 
Mod. Phys., 82, 603 

Yamada, M., Ren, Y., Ji, H., Breslau, J., Ger- 
hardt, S., Kulsrud, R., & Kuritsyn, A. 2006, 
Phys. of Plasmas, 13, 052119 

Zaqarashvili, T.V., Khodachenko, 

M. K., fc R ucker, H. O. 2011, 
http://arxiv.org/abs/1101.3913vl, 

Zweibel, E. G. 1989, ApJ, 340, 550 

Zweibel, E. G., k Yamada, M. 2009, Annu. Rev. 
Astron. Astrophys., 47, 291 



This 2-column preprint was prepared with the AAS lATJiiX 
macros v5.2. 



19 



